Entanglements in a Quiescent and Sheared Polymer Melt 
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Molecular dynamics simulations were performed for a polymer melt. In quiescent states the inter- 
chain interaction energy supported by each particle takes relatively large values persistently for long 
times if the particle is close to an entanglement. Thus, if the interaction is averaged over appropriate 
time intervals, we can detect active spots on the chains produced by entanglement constrains. If 
rapid shearing is applied, these active regions subsequently become bended in eventual zig-zag shapes 
of the chains. We also demonstrate that stress overshoot occurs with onset of disentanglement. 

PACS numbers: 83.10.Nn, 83.20.Jp, 83. 50. By 61.25.Hq 
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The dynamics of dense polymer melts has been a chal- 
lenging subject in current polymer physics |jl|,^). While 
the dynamics of short chains can be reasonably described 
by the simple Rouse model, the dynamics of very long 
chains has not yet been well understood on the micro- 
scopic level since it is governed by entanglement effects. 
The reptation theory [pip] is the most successful approach 
to date in describing the dynamics of entangled polymer 
chains. It has been supported by a number of experimen- 
tal and numerical works |p|-|To| , where experimentally ac- 
cessible quantities such as the stress relaxation function 
G(t), the incoherent dynamic scattering function, and 
the mean square displacement have been compared with 
the theoretical predictions. Understanding rheological 
properties of polymer melts is also of great importance 
and extensive simulations have been performed in this 
direction @-|ll|]. In a unique paper, Ben-Nairn et al. fl2| | 
could visualize individual entanglements using the fact 
that incidental contacts of the particles and entangle- 
ment contacts behave differently because entanglement 
constraints are long-lived. The main motivation of our 
simulation is to identify entanglements more systemat- 
ically on the microscopic level and examine how they 
influence rheology in rapid shearing. 

We used the bead-spring model B for our polymer 
melt composed of M chains with N beads. All the 
bead particles interact via the Lennard- Jones potential, 
U-Lj{r) = 4e[(er/r) 12 - (cr/r) 6 ] + e, truncated at the 
minimum distance 2 1 / f V. It serves to prevent spatial 
overlap of the particles. In addition, bonded pairs or 
consecutive beads on each chain are connected by an 
anharmonic spring potential of the form, J7fene(^) = 
~\k c R 2 ln[l - {r/R a ) 2 } with k c = 30e/er 2 and R a = 1.5a. 
As a result, the bond lengths deviate within only a few 
% from the minimum of Z/ljM + £^fene(?") given by 
0.96c. The number density and temperature were fixed 
at n = NM/V = 1/er 3 and T = e/k B - We measure space 
and time in units of a and tq — (mtr 2 /e) 1 / 2 with m be- 
ing the bead mass. We numerically solved the Newton's 
equation of motion and took data after long equilibration 
periods of order 10 6 to avoid aging (slow equilibration). 



In quiescent cases, we imposed the micro-canonical con- 
dition with the time step At — 0.005 for numerical in- 
tegration. In order to obtain accurate linear viscoelastic 
behavior of entangled polymers, very long simulations 
of order IOt^ were performed for a system composed of 
M = 10 chains each consisting of N = 250 beads, where 
Td = 6 x 10 5 is the stress relaxation time. In the presence 
of shear flow, we set At = 0.0025 and kept the tempera- 
ture at a constant using the Gaussian constraint thermo- 
stat to eliminate viscous heating. After a long equilibra- 
tion time in a quiescent state for t < 0, all the particles 
acquired a velocity jy in the x direction at t = 0, and 
then the Lee-Edwards boundary condition jR|[l4|] main- 
tained the simple shear flow for t > 0. Steady sheared 
states were realized after transient nonlinear viscoelastic 
regime. 

First, we show data of the stress relaxation function, 



G(t) = (P xy (t)P xy (0))/Vk B T. 



(1) 



at zero shear rate in Fig.l, where P xy is the xy com- 
ponent of the total stress tensor. For the shorter 
chain case of N = 25 and M = 40, G(t) is 
well fitted to the Rouse relaxation function |l5"|| , 
G R (t) = {nk B T '/N)Y£~? exp(-i/r p ), where r p = 
T i/sin 2 (pvr/2iV) with r 01 = (b 2 /12k B T 6.5 @. G(t) 
for the longer chain case of N = 250 and M = 10 
was obtained from a single extremely long simulation (~ 
5 x 10 6 = 10 9 A£). It clearly exhibits slow relaxation due 
to entanglements. Indeed, it can reasonably well be fitted 
to the stress relaxation function from the reptation the- 
ory, G rep (t) = (8nk B T/ir 2 N e )j: p=odd eM-P 2 t/T d )/p 2 , 
with N = 250, N e = 100, and r d = 6 x 10 5 for t > 2 x 10 5 . 

In order to identify and visualize entanglements for 
the case N = 250, we define the potential energy of non- 
bonded interaction on the z-th particle by 
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(*)= E UUlnW-r.m, (2) 
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where the particle j is not bonded to i. The distribution 
of Ef B (t) at each time t is nearly Gaussian, and there is 
no appreciable correlation even between adjacent beads, 
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because Ef B (t) consists mostly of rapidly varying ther- 
mal fluctuations. To eliminate the rapid components we 
examine its time average, 
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dt'E™(t + t'), 
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where the time interval r is much longer than 1 . We find 
that the variance, 
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decreases as r -1 / 2 with increasing r in the region 1 <C 
t <C Td, where (E NB ) is the average. This is be- 
cause the rapidly varying contributions behave as white- 
noises. With this time-averaging procedure, strong cor- 
relations due to entanglements become apparent in the 
non-bonded interactions along the chains. Fig. 2 displays 
instantaneous values E^ B in (a) and time-averaged val- 
ues Ef B with r = 5 x 10 3 in (b). The chains are straight- 
ened horizontally for the visualization purpose. While 
no correlations can be seen in (a), we can see active 
spots which consist of several consecutive beads having 
distinctly large values of Ef B . Then we recognize that 
there are two or three entanglements on each chain. Fur- 
thermore, to examine the correlations in Ef B along the 
chains quantitatively, we define the auto-correlation func- 
tion, 



C(n) = 1 ^ 



j-i=n 
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where the two beads i and j are separated by n on the 
same chain. The normalization factor Af is defined such 
that C(0) = 1; then, Af = MNct(t) 2 . In Fig.3, we show 
C(n) in (a) and its power spectrum, 



N-l 



P(k) = C(n)exp(-ikn), 



(6) 
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in (b) for various r in the range < t < tj. The cor- 
relation between adjacent beads C(l) is nearly zero for 
r = but increases with increasing t. The decay of 
C(n) for n < 10 becomes slower with increasing r. The 
width An determined by C(An) < 10~ 2 in Fig. 3a grows 
with increasing r roughly as An oc r a with a ~ 0.35. 
This widening arises from the motion of entanglements 
along the chains or equivalently from the motion of chains 
through tubes formed by entanglement constraints . 
At t = 5 x 10 3 , C(n) takes a negative minimum around 
n ~ 45 and a positive maximum around n ~ 80. Cor- 
respondingly, P(k) has a peak at k ~ 2-7r/80. Thus the 
average bead number N e between entanglements is about 



80, which is consistent with recent simulation results [||. 
For much larger r (> 5 x 10 4 ), however, no periodic struc- 
ture of C(n) was observed because entanglements become 
delocalized along the chains on such very long time scales. 

Next, we applied a shear flow with rate 7 = 10~ 3 (~ 
500/rd) using the same initial condition which gave the 
data shown in Fig. 2b. In Fig. 4 the chain conformations 
are projected on the xz plane at jt = 5 in (a) and jt = 10 
in (b). Because the chains are rapidly elongated, they 
eventually take zig-zag shapes bended at entanglements. 
Remarkably, the active regions detected without shear 
in Fig. 2b mostly become bended under shear in Fig. 4a. 
Here we give the same number to an active spot in Fig. 2a 
and that in Fig. 4a on the same chain if their distance 
measured along the chain remains shorter than 10. This 
coincidence unambiguously demonstrates the validity of 
the time-averaging method used in Fig. 2b to visualize en- 
tanglements in quiescent states. Furthermore, the non- 
bonded interactions in these entangled regions (mostly 
in bended parts) become increasingly amplified with in- 
creasing the strain, obviously because a considerable frac- 
tion of the stress is supported by entanglements in strong 
deformations. As a result, they can be detected even 
with much smaller r than in the quiescent condition, so 
we set t = 500 in calculating Ef B in Fig. 4. From Fig. 4a 
to Fig. 4b the movement of the entangled regions could 
be followed, so the corresponding parts in Fig. 4a and 
Fig. 4b can be marked with the same number. We can 
also see that the total number of these regions has not 
yet decreased in Fig. 4b at jt — 10, but several of them 
are approaching chain ends and arc going to disappear. 
Note that entanglements can be released only when they 
reach a chain end. Thus, the disentanglement process 
induced by shear flow is going to start in Fig. 4b. 

In the case of entangled melts, the shear stress P xy (t) 
and the normal stress P xx (t) — P yy (t) often exhibit over- 
shoot behavior in rapid shearing P, |l~o| , |l7f . In Fig. 5 
we show the shear stress P X y{t) and the normal stress 
difference P xx (t) — P yy (t) after application of shear at 
7 = 10~ 3 . They were calculated in the same run yielding 
Fig. 4. In agreement with experiments and the previous 
simulations, (i) these stress components both exhibit a 
peak around jt ~ 10 and afterwards tend to a steady 
state value and (ii) the peak of P xx (t) — P yy (t) is at- 
tained shortly after that of P xy (t). As can be seen in 
Fig. 4b, the chains are stretched without appreciable dis- 
entanglement until the stress maximum. We recognize 
that the stress components begin to decrease with onset 
of the disntanglement process. 

In summary, we have succeeded to detect and visualize 
the entanglements in a model polymer melt in both qui- 
escent and sheared conditions to obtain N e ~ 80. This 
work is supported by Grants in Aid for Scientific Re- 
search from the Ministry of Education, Science, Sports 
and Culture of Japan. Calculations have been performed 



2 



at the Human Genome Center, Institute of Medical Sci- 
ence, University of Tokyo and the Supercomputer Center, 
Institute for Solid State Physics, University of Tokyo. 
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FIG. 1. The calculated stress relaxation function G(t) for 
N — 25 and 250 (solid lines). The Rouse relaxation func- 
tion Gr{€) (doted lines) and that from the reptation theory 
Gre P {t) (dashed line) with N e — 100 are also shown. 
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FIG. 2. Distributions of non-bonded interactions on the chains with N = 250 in a quiescent state. We can see no correlations 
along the chains in snapshot values of [Ef B (t) — (i5 NB )]/cr(0) in (a). However, the time-averaged normalized interactions 
[Ef B (t) — (_E NB )]/cr(r) with r = 5 x 10 3 are distinctly strong in line segments consisting of several beads (in red or orange) 
due to entanglements. Among these segments we number those which keep to hold large non-bonded interaction values during 
rapid shearing in Fig. 4. The color map is the same as in Fig. 4. 
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FIG. 4. Snapshots of the deformed chains in the xz plane. A shear flow at 7 = 10 -3 is applied to the system with the 
configuration in Fig. 2b in the same run. The numbered segments correspond to those in Fig. 2b. Here yt — 5 in (a) and yt = 10 
in (b). The flow is in the horizontal (a;) direction, and the shear gradient is in the out-of-plane (y) direction. The non-bonded 
interactions with r = 500 are written on the chains. Use is made of the color map on the left, where the numbers are the 
normalized deviations, (Ef B - (E NB ))/a(r), with r = 500. 




FIG. 5. The stress growth functions after application of 
shear flow at 7 = 10 -3 in a model polymer melt with N = 250. 
The chain shapes and entanglements at yt = 10 are written 
in Fig. 4b. 
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